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High precision soft hadron abundance data produced in relativistic nuclear coUisions at LHC at 
y^SNN < 5500 GeV wiU become available beginning in 2007/8. We explore, within the statistical 
hadronization model, how these results can help us understand the properties of the deconfined 
quark-gluon phase at its breakup. We make assumptions about the physical properties of the fireball 
and obtain particle production predictions. Then, we develop a strategy to measure parameters of 
interest, such as strangeness occupancy 73, chemical potentials fiB and fis- 

PACS numbers: 24.10.Pa, 25.75.-q 



I. INTRODUCTORY REMARKS 

Relativistic heavy ion collisions experimental program 
has as objective the formation of the deconfined quark- 
gluon plasma (QGP) phase in the laboratory. The un- 
certainty in this experimental program is if, in the avail- 
able short collision time, 10~^^-10~^^ s, the color frozen 
nuclear phase can melt and turn into the QGP state of 
matter. There is no valid first principles answer available 
today, nor it seems, it can be expected in the foreseeable 
future. From this, and other such uncertainties about 
the QGP, arises the need to define and study its observ- 
ables, even though we are quite convinced that this is 
the state of matter that filled the Universe in its early 
stage, till hadronization occurred 10-20/iS into its evolu- 
tion. QGP is the equilibrium state in a hot Universe at 
temperature above that of lightest hadronic particle, the 
pion, kT > tUtjC? = 140 MeV (we hence use units such 
that k — c~ h~l). 

Detailed theoretical study of the properties of this new 
state of matter shows that QGP is rich in entropy and 
strangeness. These are the observables discussed here 
explicitly, and implicitly, in the context of soft hadron 
production. The enhancement of entropy S arises in the 
early stages of the collision process, because the color 
bonds are broken, and numerous gluons are formed and 
thermalized. Enhancement of strangeness s is in part 
also due to the breaking of color bonds. Furthermore, it 
is due to a modification of the kinetic strangeness forma- 
tion processes. These operate faster in deconfined phase, 
mainly because the mass threshold for strangeness exci- 
tation is considerably lower in QGP than in hadron mat- 
ter, but also because there are more channels available 
considering the color quantum numbers. 

We will not discuss the following important QGP ob- 
servables in this work, their current status at RHIC is 
discussed, e.g., in Ref. Aside from strangeness and 
entropy enhancement, another soft hadron signature is 
the shape of particle spectra, which carries information 
about the formation and evolution dynamics of the state 



of matter that is the source of these particles. Given the 
considerable increased energy, we expect a greater energy 
density in the initial stage, and thus a much more violent 
transverse outflow of matter than has been seen at RHIC. 
Such a strange transverse collective flow carries many 
particles to high transverse momenta, and produces a 
strong azimuthal asymmetry in particle spectra for fi- 
nite impact parameter reactions. Among other related 
hadronic signatures, we note a significant charm quark 
abundance, originating in primary parton reactions. The 
pattern of charm hadronization should reveal further de- 
tails about the QGP phase, just as strange hadrons do. 

Within a few years, a new energy domain will become 
accessible in study of heavy ion collisions, when the Large 
Hadron Collider (LHC) at CERN becomes operational 
in 2007. The top energy available to Pb-Pb reactions is 
Y^snn = 5500 GeV, a 27.5 fold increase compared to the 
top RHIC energy. Extrapolating the trends of SPS and 
RHIC physics, we expect a much greater entropy and 
strangeness yields at central rapidity. We will always 
address, in this work, most central head on collision re- 
actions of two A ~ 200 heavy nuclei, at LHC this will be 
Pb-Pb collision . 

At these high energies, there will be much less stopping 
power of baryon number, and thus the central rapidity 
region will be much more similar to the phase prevail- 
ing in the early Universe than this is the case at RHIC. 
One of the objectives of this work is to assess how small 
the baryochemical potential /ie can become, compared 
to /iB — 25 MeV observed at RHIC. The scale of at 
the hadronization of the early Universe is /ig ~ 1 eV 0] . 

In this paper, we address specifically the pattern of 
soft hadron production based on the assumption of a sud- 
den breakup of the deconfined hadron phase with all soft 
hadrons produced at essentially same physical conditions, 
and not subject to requirement of the absolute chemical 
equilibrium condition. We will, however, provide refer- 
ence data for the chemical equilibrium case. In the next 
section IhI we outline the statistical hadronization model 
and present its parameters. In section llTTI we discuss how 
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simple, but general, hypothesis allow to fix values of these 
parameters. We establish range of physical interest for 
strangeness phase space occupancy 7^. In section Hvl we 
develop our predictions regarding the properties of the 
hadronizing fireball. We present the range of statistical 
parameters which we expect and the physical properties 
of the fireball at its breakup. We consider particle ratios 
which could help determine the value of 7=,, which is a 
parameter in this study. In section ]V] we consider ob- 
servables sensitive to /iB, and obtain results which show 
how precise the hadron yields need to be measured in 
order to allow measurement of fiB- 



II. STATISTICAL HADRONIZATION AND 
MODEL PARAMETERS 

Statistical Hadronization Model (SHM) is, by defini- 
tion, a model of particle production in which the forma- 
tion process of each particle fully saturates (maximizes) 
the quantum mechanical probability amplitude. Particle 
yields are thus determined by the appropriate integrals 
of the accessible phase space For a system subject to 
global dynamical evolution, such as collective fiow, this 
applies within each volume element, in its local rest frame 
of reference. The SHM is consistent with the wealth of 
SPS and RHIC data available today. Systematic study 
of particle production for a wide reaction energy range 
confirms applicability of the SHM, see 

Analysis of hadron yields further facilitates a study of 
the physical properties of the hadronic fireball at the time 
of hadronization, when these particles are produced, i.e., 
undergo chemical freeze-out. A study of hadron mul- 
tiplicities, produced at energy ranging from the top of 
AGS energy to the top of RHIC energy, leads to an un- 
derstanding of the physical properties of the fireball at 
hadronization 0] . Here, we reverse the approach — using 
the systematics of the energy dependence of the physi- 
cal properties of the fireball, we establish our expecta- 
tions about the statistical parameters and relative parti- 
cle multiplicities expected at LHC. We cannot address, 
in this work, the total hadron yield, since this depends 
on the early stage of the reaction, and specifically, on the 
entropy formation process in initial interactions. 

The complexity of the SHM model derives from the 
need to account for many hadronic resonances, and their 
decays. Since the number of massive resonances grows 
exponentially, their contributions to particle yields, es- 
pecially to the yields of pions, are slowly convergent. 
Moreover, the counting of resonances and their conserved 
quantum numbers poses a significant book keeping chal- 
lenge. For this reason, an effort has been made to gen- 
erate a comprehensive software package available to all 
interested parties, with a transparent hadron data in- 
put, and a comprehensive parameter field. All results 
presented were obtained using this numerical package 
SHARE (Statistical Hadronization with REsonances) [a| • 

Another package of similar capability, but with re- 



strained parameter set, has since become available 

In- 

These programs are including a large number of reso- 
nances and track the chemical composition as well as the 
decay trees with care. As result, the benchmark fits pro- 
duce a hadronization temperature which is considerably 
lower than obtained in Ref. 01 for SPS or RHIC reaction 
systems. 

A successful description of rapidity particle yields 
within the SHM, at a single-chemical freeze-out condi- 
tion, produces the model parameters in the process of 
minimization. The parameters are: 

1) dV/dy, the volume related a given rapidity to the par- 
ticle yields; 

2) T, the (chemical) freeze-out temperature; 

3) = Tln(A„A(j)'^/^, the baryon and 

4) lis = Thi[Xq/Xs], hyperon chemical potentials; 

5) A/3 = Xu/Xd, a fugacity distinguishing the up from 
the down quark flavor; 

6) 7s the strangeness phase space occupancy; 

7) 7g the light quark phase space occupancy. 

When the assumption of absolute chemical equilibrium 
is made the values 7^ = 7^ = 1 are set. The relationship 
to quark flavor fugacities Xi,i = u,d,s is made explicit 
above. 

When a set of parameters is known, all particle mul- 
tiplicities can be evaluated exactly. Similarly, one can 
obtain the physical properties of the fireball such as ther- 
mal energy, entropy, baryon content by appropriate eval- 
uation of the properties of partial fraction contributions 
of each hadronic state. In that way, in fact, one obtains 
from a fit to a limited set of measured particle yields 
a full phase space extrapolation for any particle yield, 
and a full understanding of the properties of the fireball. 
Some striking 'conservation of fireball properties' rules 
emerged from such an analysis of particle yield data 0| , 
and these we will use in order to predict the (relative) 
particle yields at LHC energy range. 

III. LHC HADRONIZATION 
A. Choice of conditions and constraints 

To predict the 7 parameters, at first sight, we need at 
least 7, and better more than that, valid conditions, con- 
straints and hypothesis: 

1) The 'volume' normalization dV/dy only enters abso- 
lute hadron yields. Consequently, it is related to the ini- 
tial conditions, i.e., mechanisms of entropy production. 
Restricting our investigation to the study of particle ra- 
tios, we do not need to know the value of this parameter, 
which normalizes the overall yield. 

There are two natural physics constraints: 

2) Strangeness conservation, i.e., the (grand canonical) 
count of s quarks in hadrons equals s count at each ra- 
pidity unit. In our specific case, we request that: 

^-^ = 0±0.01. (1) 

s + s 



3 



We will show how this condition establishes the relation- 
ship between the baryochemical potential /zb and the 
strangeness chemical potential fis in some detail in sub- 
section 

3) The electrical charge to net baryon ratio, in the final 
state, should be the same as in the initial state, and in 
the specific case of Pb-Pb interactions, we have: 

y = 0.39 ±0.01. (2) 

The 2.5% error can be seen as expressing uncertainty 
about how well neutron and proton densities follow each 
other, given that even most central reactions occur at fi- 
nite impact parameter, and 10% of nucleons do not par- 
ticipate in the reaction. 

4) and 5) 

i) The phase space occupancy parameters in some ap- 
proaches are tacitly fixed: assuming the chemical equi- 
librium one sets 

7s = 79 = 1 • 

ii) In the chemical non-equilibrium approach, the sys- 
tematics of data analysis at RHIC and high SPS energies 
firmly predicts: 

Furthermore, we will present our results as function of 
7s. We believe that the value of 7^ is linked to the colli- 
sion energy, as more strangeness can be produced, when 
the initial conditions reached become more extreme. We 
show, in section HVBI how 7^ fixes several easily accessi- 
ble observables and can be measured, and the consistency 
of the chemical non-equilibrium approach checked. At 
LHC energy, the expected value of 7^ is so much greater 
than unity, and thus, we can be sure that a distinction 
from 7s = 1 for the equilibrium model can be arrived at 
in an unambiguous way. 

There are (at least) two further conditions required 
which must be sensitive to the hadronization temperature 
and baryochemical potential. These will be drawn from 
the following observations: 
6) We note the value, 

^-.0.78, or ^ 0.845, 

for chemical non-equilibrium Q, or respectively, equi- 
librium model analysis. We note that the energy per 
particle of non-relativistic and semi-relativistic classical 
particle gas is E/N ~ m -f 3/2 T -I- . . ., while the entropy 
per particle in this condition is S/N = 5/2 + m/T + . . . 
(see section 10 of 0). Hence: 

E_ ^ m/r + 3/2 

TS m/T + 5/2 ■ ^ ' 

It is thus possible to interpret this constraint in terms 
of a quark matter made of particles with thermal mass 



TO cx aT. Solving Eq.El we find for E/TS = 0.78, a = 2 
for chemical non-equilibrium, and for E/TS = 0.845, 
a — A for equilibrium. This is near to the result ex- 
pected in finite temperature QCD |[8|. That result points 
to a simple structure of the quark matter fireball. On 
the quark-side, the value E/TS is not very model spe- 
cific, though it is sensitive to the average particle mass 
as shown above. 

On the other hand, there is considerable sensitivity to 
this thermodynamic constraint in the hadronic gas. The 
hadron system after hadronization comprises a mix of 
particles of different, and for the baryon component, large 
mass. To fine tune this value a specific ratio of baryons to 
mesons needs to be established: in this way the hadron 
system can maintain both energy and entropy, aside of 
baryon number and strangeness, during hadronization. 
For this reason, in the hadron phase there is consider- 
able sensitivity of this ratio to both T and the value of 
the phase space occupancy parameters, here 7s. 
7) We need to make an assumption which fixes the bary- 
ochemical potential /ib- This certainly is one of most 
difficult guesses as there is no reliable way to predict 
baryon stopping at LHC, and certainly this value, 

will be quite difficult to measure. For this reason, we will 
discuss, in subsection FVl at length a method to measure 
/.tB. The impact of /iB — 1~3 MeV on mixed particle 
ratios, such as K/tt or jy/Tr" is physically irrelevant. We 
will show haw particle-antiparticle abundances can help 
us further. In fact, it is uncertain that a measurement 
of such small /zb can be accomplished, and thus in ef- 
fect, we could have simply assumed /zb = which at a 
few %-level would be consistent with all relative hadron 
yield predictions here presented. On the other hand, the 
understanding of matter-antimatter asymmetry present 
at LHC energy scale is, in itself, of interest and thus, we 
pursue this question further. 

We argue as follows: our analysis of RHIC data sug- 
gested that baryons are more easily retained in the cen- 
tral rapidity region than energy, with 2.5 times larger 
fraction of colliding baryons than fraction of energy de- 
posited. Considering the SPS data point, and the RHIC 
results, the per baryon thermalized reaction energy re- 
tained in the central rapidity region, in units of the max- 
imum available collision energy, drops from 40% avail- 
able at RHIC to 15% at LHC. Given the LHC energy 
■y/SNN = 5500 GeV, a thermal energy content per net 
baryon at y = is assumed to be, 

dE 

— = 0.15 X 5500/2 = 412 ± 20 GeV, 
do 

at top LHC energy, where the error is chosen to be similar 
in relative magnitude as the error in other observables 
considered. This error plays a role in finding the solution 
in terms of statistical variables of the constraints and 
conditions considered. 
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This specific assumption fixes implicitly the (small) 
value of the baryochemical potential, and by virtue of 
the strangeness conservation, also of the strange chemi- 
cal potential. When we deviate from this assumption in 
exploring a wider parameter range, we will mention this 
explicitly. 



B. rj and maximum value of 73 

We will study the hadronization condition as function 
of 7<; which may take large values. We note that 7^ can- 
not rise above a limit, to be determined from similar 
consideration as is the maximum value of 
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exp(TO7r/2r). 



(4) 



At this value, the Bose distributions of pions diverges. 
As 7s increases, same will happen in the strange hadron 
sector and indeed this will occur first to the lightest 
particle with considerable hidden strangeness content, 
i.e., ?7(548). Naively, one could expect that 7^^ = 
exp(mj,/2T). However, 77, unlike the spin 1 0(1020) is 
not a (nearly) pure ss state. 

The quark structure of ry(548) and 77'(958) can be writ- 
ten as: 



dd 



^/2 _ 
uu + dd , 



(5) 



^/2 



Study of numerous experimental results, and in particu- 
lar of Z'^ — ^ hadrons, shows that sin^cpp — 0.45 ± 6 indi- 
cating that 77(543) has 45% strangess content Q. This 
arises from the SU(3)-flavor octet state content of 67% 
reduced by SU(3)-symmetry breaking mixing with the 
33% strangeness content of the singlet 77' (958). 

In order to count the yield of the 77, we introduce its 
fugacity T^. The fractional contribution to the partition 
function is: 



In Zri 



dV- 



d^p 



(27r)3 



In 1 



(6) 



with E,i = y 777^ + p2 focus our attention on the 

dominant, directly produced 77. The incremental per unit 
volume rj yield is: 



dV " " 



d[d\nZ^/dV] _ f d^p 1 



9T„ 



(2^)3 T-ie*-l' 



More specifically, the rapidity density is: 

dN.„ dV dN„ 

rir, = — X -. 

^ dy dy dV 



(7) 



(8) 



The 'normalization' dV/dy which comprises the trans- 
verse dimension at hadronization, and the longitudinal 



incremental volume, is arising from kinetic expansion 
processes driven by the initial state formation mecha- 
nisms. These have to be obtained by methods beyond 
the scope of this work. However, we note that in the lon- 
gitudinally scaling limit for ideal fluid hydrodynamic evo- 
lution of the initial state, at all times the entropy rapid- 
ity density remains constant, dS/dy — Const.. Since the 
particle multiplicity is defined by the value of dS/dy, the 
total hadron multiplicity is not affected by our ensuing 
study of hadronization for different chemical freeze-out 
temperatures. The total charge particle rapidity density 
is a consequence of initial processes which are beyond the 
physics reach of this work. 

We now relate the ?7-fugacity T,, to the light and 
strange quark fugacities 7, and 7^. In the SHM, the 
probability of the production of 77 is weighted with the 
yield of strange and light quark pairs in proportion of 
their contribution to the quark content in the particle 
formed. Thus, 



T^=72cos2 



2 ■ 2 

-7sSm I 



< e' 



(9) 



The upper limit is set at the phase space divergence 
point. Specifically, considering 77(548) with its 45% 
strangeness content, we find from Eq. a maximal value 
7s < 10.4 for the hadronization conditions of interest in 
this work, i.e., 7, e™''/^^ and T 140 MeV. Thus 
we set as the range of interest < 7s < 10. 

This upper limit is the most stringent constraint for 
7s: it is more stringent than the one which arises from 
the consideration of the (A-meson, 



(10) 



due to the greater (jj-mass. Similarly, the constraint on 7s 
based on the kaon condensation (in presence of negligible 
chemical potentials): 



Tk = 7s 7<? < e 



(11) 



is less severe. 

To confirm the functional dependence on 7^, Eq. (|5Jl, 
we show that it is consistent with the requirement that 
the fugacity 7^ , allows the count of valance quark content 
in hadrons. We look at 7s, which allows the count of all 
strange and antistrange quarks by the relation: 

dlnZ 

s + s = 7s— . (12 

The contribution of 77 to the strangeness count thus is, 

(, + ^\ - 2±^r ^^^^^ 

^'^''^ - T, 97. " dT, ' 
o 7sSm^0P 



72 cos2 (Pp + 7^ 

which is the required result. A similar procedure can be 
followed to show that: 



2 Qin^ , 



-N^, (13) 



cos^ 0p 
7g cos2 (\)p + 72 sin^ 



(14) 
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We conclude that Eq. ^ defines the 7/-fugacity in terms 
of light and strange quark fugacities. 

C. Range of 73 of interest 

When 7s increases, the 77-fugacity increases rapidly, 
but, even at 7s = 5, it is well below the Bose singularity. 
For T — 160 MeV, the singularity would be just above 
7s = 8. However, the fall-off of the expected hadroniza- 
tion temperature (see below) moves this towards twice 
as large value. One may wonder how large a value of 7s 
is physically consistent, and in particular, if an associate 
decrease in the value of T makes sense. 

In qualitative terms, this type of parameter evolution 
and correlation is fully consistent with the picture of 
rapid transverse expansion of the QGP fireball. This 
expansion can lead to supercooling which pushes the 
hadronization temperature lower. At the same time, the 
preserved yield of strangeness requires that 7s increases. 

We now look at the possible range of expected val- 
ues of 7s in more detail. It can be expected that, for 
LHC extreme conditions, the strangeness phase space has 
been chemically saturated at a temperature larger than 
strange quark mass: Ti > rUs- Moreover, in the QGP 
phase, there is a residual strange quark mass mg > T, 
where T is the final hadronization temperature. While 
the precise value of ms depends on the momentum scale 
at which it is measured, we will assume in this semi- 
quantitative discussion that ms(l GeV) ~ 180 MeV. 
Conservation of strangeness rapidity yield in the expan- 
sion from Ti to T implies, by comparison of the relativis- 
tic phase space size, 

{dViTi)WimjTi) = 7sQ X {dVT^)Wims/T), (15) 

where 7^ is the quark-side hadronization phase space 
occupancy, for definition of W{x), see Eg. (|18f) . For 
< X < 1, W{x) is slowly changing and close to its 
asymptotic value W{x = 0) — 2. Entropy conservation 
further requires that dViT^ = dVT^ = Const. 

For Ti ~ 1.5r, and noting that Ti > nig > T, we ob- 
tain, as solution of Eq. (|15() . 7^ ~ 1.5. In the hadroniza- 
tion process, this value increases by a significant fac- 
tor which we can obtain comparing the phase space of 
strange hadrons in QGP with that hadron phases (we 
omit as is customary the upper index H, i.e., 7^ — > 7s): 

-lfxW{m,/T) = 7sX^VK(m,/r) (16) 

i 

+ 7s x5]VF(m,/r) 

3 

+ 7s'x5]w^(™fc/r). 

k 

The sums run over single, double and triply strange (s 
and s) hadrons, respectively. W differs from W in that it 
comprises appropriate hadron fugacities for each particle. 



the appropriate expressions are seen in more detail in 
Eq. (|T7|l. 

The strangeness QGP to HG aspect ratio, i.e., the ra- 
tio of the phase phase space size seen in Eq. (|16|l up to 
the coefficients 7P, 7s shows that it grows with decreas- 
ing temperature, see Fig. 19.3 in Ref. This is further 
strongly amplifying by a factor as large as four the final 
observed 7s. A QGP phase value 7^ ~ 1.5 may be- 
come 7s — 5-10 on hadron side. As we shall see, the 
decrease of T with increasing 7s compensates, in part, 
the increase with 7s in the yields of strange hadrons. We 
will thus consider the range 0.5 < 7s < 10 as that is 
where most of variation in the considered observables is 
seen, and the actual physical conditions are expected to 
occur. The lower limit underscores the comparison with 
the equilibrium model 7s = 7g = 1 behavior. The upper 
limit 7s < 10 is within the range of allowed values of 7s 
we have obtained in previouse section IlII Bl 



IV. PREDICTIONS 
A. SHM parameter values at LHC top energy 

With the assumptions outlined in subsection IIII Al we 
solve for the best set of SHM model parameters. A pre- 
cise solution is always found in the considered range of 
7s, also for the increased value of E/TS for the case of 
chemical equilibrium. This of course does not guarantee 
that, e.g., our baryochemical potential is correctly cho- 
sen. The results we show are not entirely smooth as we 
allow a small error in the constraints and conditions, and 
thus the solution for the parameters is not a precise alge- 
braic result, but a most likely value of a quasi-fit, which 
has a ~ 0. 

The resulting statistical parameters T, 7^^, /ib, and 
/j,s are shown, in Fig.^ as function of 0.5 < 7s < 10. 
Note that through our study of chemical non-equilibrium 
7q ~ 7^^. We see that an increasing value of 7s is accom- 
panied by considerable reduction of the hadronization 
temperature, which drops from the value T = 140 MeV, 
near 7s = 2.4 and 7, = 1.6, down to T = 110 MeV. At 
the favorite value 7s = 5, the expected LHC hadroniza- 
tion temperature is T = 125 MeV for the chemical non- 
equilibrium. The chemical equilibrium result at T = 156 
implies a chemical potential /ib = 2.6 MeV, and this 
value is also found for the chemical non-equilibrium 
model for 7s = 5. We explain why /xb remains un- 
changed at the end of section For , we see a slight 
reduction by 10% from equilibrium model value at 0.52 
MeV as is appropriate considering the results presented 
in Fig.[Sl Our estimate of the chemical parameters at 
LHC are considerably different from those proposed by 
others where ^b = 1 MeV and /xs = 0.3 MeV is 
proposed. We note that our baryochemical potential is 
considerably greater. We will discuss experimental con- 
sequences at greater length in section IVl 

It is of considerable interest to study the physical 
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FIG. 1; (color online) The values of T, 7^ , /^b, and fis 
as function of varying 7^, consistent with the hadronization 
model assumptions outlined in subsection IIII Al The equilib- 
rium model results are crosses at 7^ = 1 for 7, = 1. 
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FIG. 2: (color online) Pressure P [GeV/fm'^], energy density 
e [GeV/fm'^], entropy density a = S/V [1/fm^], net baryon 
density f = {B -B)/V = b/V [1/fm^], for non-equilibrium 
SHM. Cross at 7^ for chemical equilibrium. 



properties at hadronization: the pressure P, the energy 
density e — E/V, entropy density a = S/V, and net 
baryon density ly = B/V. In general, the non-equilibrium 
hadronization occurs from a state of higher density, as is 
seen in Fig.|21 comparing the lines with the SHM equi- 
librium cross. This is, in particular, true for the entropy 
density. With increasing 75, all density decreases, the 
drop in temperature is a more important influence than 
the large increase in relative strangeness yield. Despite 
a significant increase in fi^ with increasing 7^, the net 
baryon density decreases modestly. It is very small, bor- 
dering the value 1/ = 0.001 fm~'^. 



B. Particle yield ratios and determination of 7^ 

There is approximate charge symmetry with positives 
and negatives h~ having a very similar yield. The dif- 
ference will be discussed further below. The total charged 
hadron yield will be denoted as, 

h = h'^ + =p + p + 7r+ +71^ +K+ +K~, 



and is evaluated after weak decay of hyperons and Ks,l- 
Similarly the total yield of neutrals: 

h° = 11° + n + n. 

The top panel, in Fig. 13 shows the 2hP/h ratio, which 
varies by ±10% in the range of 7^ considered, with the 
charge symmetric value 2hP /h = 1 arising at 7^ = 3.3. 

In the lower three panels, we focus our interest on ra- 
tios of some interest for determination of 7^. We present 
the ratio rj/i:^ which rises by nearly 50% compared to 
the equilibrium model expectation, see crosses at 7^ = 1, 
computed for 7g = 1 and T = 156 MeV. This ratio is ob- 
served by reconstruction of the invariant di-photon mass. 
This observable shows a relatively small variation with 
7s , which can be understood as result of competition be- 
tween 7s and 7^, see Eq. which is accompanied by the 
decrease of T with increasing 7s . We record that the ex- 
pectation for the non-equilibrium hadronization at LHC 
is: 

0.07 < 4t < 0.12. 

This interesting observable may also not have the sensi- 
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FIG. 3: (color online) The predicted yield ratios as function 
of 7s, from top to bottom: ratio of neutral to charged hadrons 
2h°/h, r)/-K°, (K+ + K-)/h and 24,/h as function of 7^. The 
cross indicates chemical equilibrium model prediction. All 
yields after weak decay of hyperons and Ks,l. 



tivity required to distinguish the models, or help deter- 
mine 7s. 

In the panel below, we show the (K+ + K^)/h ratio. 
This ratio is rising, for large 7^, to a value near 0.23 
almost double the 'standard' chemical equilibrium value 
at 0.12. However, for 7^ = 5, we see a more modest 
increase by 50%. We record for LHC: 

K+ +K- 
0.12 < , ^ , < 0.23. 

h+ + h- 

More spectacular is the expected increases in the 2(j)/h 
ratio. The Chemical equilibrium value of 0.015 is seen 
to rise 4-fold, and at 7s = 5, we still see a very notice- 
able increase by a factor 2.5. This is a very important 
observable of the condition of hadronization. We record 
our LHC expectation: 



0.015 < 



20 



< 0.06. 



The behavior of the baryon yields is shown in Fig.0] 
where the total, nearly matter-antimatter symmetric 
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FIG. 4: (color online) (Strange) baryon ratios with charge 
hadron multiplicity (p + p)/h, A + A)//i, (H^ + H )/h and 
(yi~ -\- Q, )/h, after all weak decays occurred. Crosses denote 
chemical equilibrium result. 



yields of protons, and the three hyperon families are 
shown, normalized by the total charged hadron yield h. 
The lines, from top to bottom are for {p+p)/h, {A+A)/h, 

(S~ +E )/h and (fi" +^ )/h. We note that for large 
7s, a considerable change in the expected baryon popu- 
lations ensues, with proton yield decreasing (due to de- 
crease in T while the more strange the hyperon is, the 
more it is enhanced compared to chemical equilibrium 
expectations. Interestingly, these relative yields saturate 
with increasing 7s, as the effect of temperature decrease 
competes with the increase due to rising 7s. Note that 
we have shown, in Fig.^ the total hadron yields after 
all weak decays have occurred. We summarize our LHC 
expectations: 



0.07 > 



P + P 

h+ + h- 



> 0.04, 



0.02 < , A + f < 0.04, 



0.004 < 7 , + - < 0.015, 



0.0006 < 



h+ + h- 



< 0.004. 
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V. MEASUREMENT OF THE 
BARYO-CHEMICAL POTENTIAL 

A. Strangeness conservation 

An interesting challenge, at LHC, will be the mea- 
surement of chemical potentials. We recall that Ms is 
directly related to fj.^ should hadron emission at each ra- 
pidity occurs such that there is local strangeness balance, 
i.e., 'conservation'. This relates the two chemical poten- 
tials as we shall next discuss. Otherwise, if emission of 
hadrons were to reflect on the QGP conditions, we would 
expect /zs = and this would generate a buildup of resid- 
ual strangeness in a distillation process 113. We tacitly 
assumed, and continue in this way now, that at LHC, 
in each rapidity region, local conservation of strangeness 
prevails. We present a set of results which will allow, 
given appropriate experimental sensitivity, to determine 
the value of Ms- 

Strangeness conservation establishes a relation be- 
tween the chemical potentials While this relation- 
ship is simplified for the case ^J.^/T 1, the presence of 
7s 3> 1 introduces new elements and we reinspect the re- 
lationship. Using, at first, the quark fugacity notation for 
convenience, the open strangeness sector partition func- 
tion is, in the here appropriate Boltzmann approxima- 
tion: 



\nZs = "/sIqFk 



A2A, 



itFn ( + - ) . (17) 



The phase space factors, 



VT' 



2tt 



comprise all contributing hadron states 'fc' with quantum 
number H\ W{x) — x^K2{x) 2 is the relativistic 
phase space integral in classical limit, which for large 
X behaves as W{x) cx x^/^ exp(x/r). A series of these 
terms represents quantum statistics, see Eq. (|19|l . 
The strangeness conservation condition. 



X.dlnZ 



= 0, 



yields for small values of chemical potentials, in units 
of T: 

MB- MS + 2(mb- 2ms)-|^ + 3(mb- 3Ms)4f^ (19) 



Fk 



MS 



7sFK(2mK) 7<j7s^K(3mK) 



IqFy 



9Fy 



The right hand side presents the three first terms of the 
kaon Bose integral expansion which keeps terms of the 



20 



10 



CO 





\\ ^ 






\\\ ^ 












VsX \ 






\ \\\ ^ 






\ \V\ ^ 






\ \\\ 






\ \ \x\ ^ 






\ \ \\\ ^ 






\ \ \ ^ 






\ \ \ \x \ ^ 






\ \ \ \ \\.^\ ^ 
























. . 1 . . . 1 . , . 1 


1 1 1 



100 120 140 160 

T [MeV] 



180 



FIG. 5: (color online) /xb/ms ss function of T: From top 
right to left 7^ = 0, 0.5, 1, 3, 5, 7, 10, at 7q — exp(m^o/T). 
The dashed (red) line shows the chemical equilibrium model 
result at 7a = 1 and 7, = 1. Crosses correspond to RHIC 
freeze-out conditions both in equilibrium (at dashed line) and 
non-equilibrium. 

same order as those belonging to S and ^ in the balance. 
These are not entirely negligible for large 7s. We solve 
Eq. ({T^ and show, in Fig.|Sl mb/ms as function of T: 



— = fiT;7q,is)- 
Ms 



(20) 



The dashed curve is the equilibrium case for 7s = 7q = 
1, where we marked the RHIC hadronization condition 
Mb/ms = 5 with a cross. The closest solid line to this 
result is for 7^ — e™""/^ and 7s 0. The following 
solid lines are, in sequence from upper right, for 7s — 
0.5, I, 3, 5, 7, 10. The cross near the 7s = 3 line, at 
Mb/ms — 4.6, corresponds to the RHIC chemical non- 
equilibrium hadronization condition. 

These results allow, aside of gauging the LHC values 
of chemical potentials, also an easy comparison with and 
cross check of other work addressing LHC and RHIC en- 
vironments with conserved strangeness. We recall that 
many particle ratios are directly determined by chemical 
potentials, e.g., K+/iC~ = exp(2/ig/r) and S^/S~ = 
exp((2/LiB — 4:fis)/T). This obviously leads to consistency 
conditions in the particle-antiparticlc asymmetry sensi- 
tive particle ratios. 



B. Hadron— ant ihadron asymmetry 

In order to measure chemical potentials, we need to 
be able to measure the particle-antiparticlc asymmetry. 
This is not an easy task at LHC as we shall see. Up to 
RHIC energy, it was customary to study ratio of antipar- 
ticle yields to particle yields, such as A/A. At LHC, the 
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FIG. 6: A/i as function of 7s. Cross indicates equilibrium 
model result. 



strategy has to slightly change. We consider of normal- 
ized particle-antiparticle difference yields: 



1 1 / 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



A7V,: = 



(21) 



where Ui is the rapidity density dN/ dy of charged parti- 
cles i — K+, p, A, 5, and antiparticles. We omitted 
intentionally the pion from this list, as it is the dominant 
component of the unidentified particle hadron asymme- 
try, we will address next. Below, we will omit the factor 
d/ dy as the expressions we state are valid more generally. 
The most accessible observable. 



A/i = 



h+ - h- 



(22) 



composed of unidentified charged hadrons is very hard to 
measure precisely, for i) there are distortion possible by 
partial acceptance of weak decays, and ii) this variable 
assumes a comparatively small value at LHC. Moreover, 
this is also an observable hardest to interpret in a sim- 
ple and transparent theoretical model. We find that its 
magnitude for LHC will be: 

A/i = (0.5-0.6)10"^ 

for all values of 7^ > 1, as is shown in Fig.|Sl The dif- 
ference in A/i between equilibrium and non-equilibrium 
result is due to the increased entropy and thus hadron 
multiplicity content of the chemical non-equilibrium case. 
The variable A/i has been recognized already in the study 
of SPS reactions as a sensitive probe of entropy produc- 
tion [Hill. 

Should a way to measure A/i of this magnitude with 
reasonable precision be found at LHC, this would pro- 
vide a model independent measure of the value of the 
baryo-chemical potential /ib- To see this, we solve (fit) 
the SHM with a fixed assumed value of A/i. In the chem- 
ical equilibrium version of SHM, we find, at 7s = 7q = 1, 
in the T-ji-Q plane the constraint lines determined by as- 
sumed values of A/i, shown in the top panel of Fig.[7| 
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FIG. 7: (color online) A/i constraint from left to right for 
LHC Ah = 0.001; for RHIC = 0.01 and for SPS top energy 
= 0.1. Top panel: chemical equilibrium SHM, T-fiB plane, 
bottom three panels: chemical non-equilibrium, from top to 
bottom T, followed by 73 and dE/db combined with hb- In 
the bottom panel the dashed line indicate the systematics 
of the behavior regarding the value of thermal energy per 
baryon. Vertical dotted lines brace the extreme allowable 
values of fiB- 



From left to right, we have considered three benchmark 
values. Ah = 0.001 (top LHC energy), Ah = 0.01 (top 
RHIC energy), and Ah = 0.1 (top SPS energy, central 
rapidity region). At SPS, top energy in S-Pb interac- 
tions, we have studied this variable (called Dq ~ 0.009) 
in the context of the exploration of the specific entropy 
per baryon |l2l Il3j | , and the results we present here are 
in agreement with this early study of SPS hadron multi- 
plicities. 

For non-equilibrium case, with 7^ = 7^^, the corre- 
sponding result is shown in the second top panel. The 
nearly vertical lines are solution of the the following 
conditions imposed: strangeness conservation, charge-to- 
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baryon ratio and Aft, as function of /ib- For each value 
of T, there is a specific value of 7^ indicated in the next 
lower panel, and again, we see nearly a vertical line. Both 
these results imply that while the value of T and 73 are 
not much constrained by the measurement of just one 
single observable Ah, the value of /ib is already highly 
constrained. 

Why this is the case is understood inspecting the bot- 
tom panel in Fig.[7| We show the resulting value of 
thermal energy per baryon E/b. These turn out to be 
highly localized regions. Thus, Ah is for a wide range of 
other statistical parameters closely related to the value of 
energy per baryon, or equivalently, entropy per baryon. 
The connecting dashed line in the bottom panel guides 
the eye. A check of E/TS also confirms that these so- 
lutions produce the expected result which is otherwise 
introduced as a constraint. 

Thus, we learn that just the single 'measured' value 
Ah is enough to constrain a rather narrow range of /ie- 
This result is indicated by the vertical dotted line, which 
we place bracing the domains of that are allowed. 
Inspecting the equilibrium model intercept, we realize 
that this singles out a domain of T which is result of 
data fits with this constraint. This explains why /iB is 
usually determined in a model independent way within 
the SHM, with little if any difference present between the 
different model variants, provided that the experimental 
data used in the fit comprises explicitly, or implicitly. 
Ah. For example, the study of the impact parameter 
dependence at RHIC using different SHM model variants 
produced /ib and which cannot be distinguished (see 
bottom panel of Fig. 1 in Ref. 0]). 



C. Identified particle— antiparticle asymmetries 

For the identified particles, the normalized particle- 
antiparticle differences can be closely and analytically re- 
lated to the value of chemical potential. For example, the 
kaon asymmetry is directly related to strangeness chem- 
ical potential /is: 



AK 



K^ 



K^ 



K- 



tanh 



MS 
T 



T ■ 



(23) 



We have, for simplicity, not considered the (/)-meson decay 
contributions which increase the normalizing yield but do 
not alter the difference. 

We show the actual, with all decays, AK as thick solid 
line in the top panel of Fig.|Sl (bottom line in this top 
panel). The very top short dashed line in the panel (red 
on-line) is the ratio fJ,s/T. The thick long-dashed line 
excludes from the ratio the contamination by the decay 
(j) —^ K"*" -|- K~. The fully weak decay contamination cor- 
rected results are the thin (solid and resp. dashed) lines 
at the top of the parallel lines, and are shown for both 
the full result (solid lines) and </)-decay corrected result 
(long dashed). The parallel line regions are where the ac- 
ceptance of weak decays is partial and/or the correction 
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FIG. 8: (color online) Solid lines: the relative particle an- 
tiparticle asymmetry as function of 7^, from top to bottom 
AK, Ap, AA and together in bottom panel AH and Af2. The 
short dashed line (red) is the analytical result, given in terms 
of i-ii/T (see text). The range of weak decay corrections is 
shown by parallel lines with the most decay contaminated re- 
sult being the bottom, thick line. For AK we also show by 
a long dashed lines the result after the 0-decay contribution 
to AK has been removed. The equilibrium model results are 
shown as crosses (after weak decays) and circles (corrected for 
weak decays) at the 7s = 1 vertical dotted line. The diamond 
symbol in bottom panel is the chemical equilibrium result for 
AH, different from the cross for AO. 



is incomplete. After the removal of the (jj-decay dilution 
of the kaon yields, Eq. (|23ll should read: 



AK ~ 0.9—. 

T 



(24) 



The slight reduction from the analytical formula Eq. 123|) 
is due to the strong decay contributions of hyperon res- 
onances decaying emitting a kaon. Because of the small- 
ness of AK, the baryon asymmetry in the hyperon reso- 
nances leaves this small but visible imprint of this result. 
The chemical equilibrium result (cross before and cir- 
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cle after weak decays) is also indicated in Fig.|Sl These 
are, in general, larger than the analytical results (short 
dashed lines) except in the case of Aif . 

For baryons there are four particle-antiparticle differ- 
entials, which are shown below AK, in Fig.|S| We expect 
for protons: 



Ap^^^tanh^^^. 
p + p T T 



(25) 



The thin solid line in the second panel from the top which 
corresponds to removed weak decays compares well to the 
analytical results, short-dashed line. The thick solid line 
at bottom of parallel lines includes in Ap the contami- 
nation from weak decays of A and A, and the region in 
between spans all possible WD contamination. 

Once the weak decay contributions 51, S — > A and 
17, S ^ A are removed, we have further: 



AA 



A 



A + A 



= tanh 



T 



Mb - Ms 
T 



(26) 



The thin line, in the AA panel, is nearly indistinguishable 
from this result. The thick solid line includes all weak 
decays. 

There is no significant contamination of and S 
and thus we have: 



AS = 



= tanh 



and similarly for the il 



An 



n- + n 



= tanh 



T ^ T ' 



Mb - 3^s Mb - S/xs 
T ^ T ' 



(27) 



(28) 



both shown in bottom panel of Fig.|Sl 

For baryons, the expected asymmetry is at %-level and 
the weak decay of hyperons allow unique identification of 
these particles. It is quite possible that the measurement 
of these variables will succeed. As the above expressions 
show, there is a consistency condition since two chemi- 
cal potentials and the temperature considering that the 
isospin bath of numerous pions causes A/3 ^ 1 to a great 
precision, such that Tin A/3 ^ us determine five ob- 
servables AA''. This consistency is further tightened due 
to strangeness conservation relation of mb/ms = 4-5, see 
section lV Al The strangeness conservation constraint and 
through it the SHM model can be tested by transforming 
Eqs. H23I25I2()I27I28|) so we can make use of Eq. lO and 
the values Eq. (pn|l : 



Ap 
AK 
AA 
AK 
AS 
AK 

An 

AK 



Mb 
MS 

Ms 



MB _ ^ 



^-2 
Ms 

Mb 

Ms 



4.5, 

~ 3.5, 
~ 2.5, 
3 ] ~ 1.5. 



(29) 
(30) 

(31) 
(32) 



These relation test the SHM and strangeness conserva- 
tion, they do not differentiate model variants such as 
chemical equilibrium and chemical non-equilibrium. 



VI. FINAL REMARKS 

The total multiplicity yield, as well as the yield of 
charmed particles, is originating predominantly in the 
early stage, primary parton reactions. For this reason, 
we did not address absolute yields of hadrons, and, sim- 
ilarly, cannot study the total charm yield in the context 
considered here. However, one may wonder, if the ap- 
pearance of small but distinguished charmed meson and 
baryon yield does not offer an interesting and indepen- 
dent probe of the properties of the hadronization state, 
or even it could influence the results presented. 

A earlier study of charmed hadron production, in 
chemical equilibrium at T = 170 MeV, has yielded in- 
teresting insights into the relative production strengths 
of charmed mesons and baryons emerging at this partic- 
ular hadronization temperature |lOj| . We will not enter 
into further discussion of this subject, but note that: 

a) for 7s > 2, these results imply that Ds(cs)+ and its an- 
tiparticle should be the dominant charmed hadron frac- 
tion; 

b) the differences in yields between particles and antipar- 
ticles, seen in table 3 of Ref. jlfl], indicate that charm 
particle contribution to the asymmetries we studied are 
totally negligible. It appears that measurement of rela- 
tive yield of charmed mesons and baryons will reveal the 
charm hadronization condition, and we hope to return to 
this subject soon. 

The small relative number of charmed quarks and 
their even smaller particle-anti-particle asymmetry as- 
sures that these particles do not impact any of the re- 
sults we obtained. On the other hand, these particles 
offer another opportunity to explore hadronization con- 
ditions. The formation of charmed hadrons is expected 
to occur prior to general hadronization, considering the 
greater binding of charmed particles |l5| . 

To summarize, we have presented a detailed study of 
the soft hadron production pattern at LHC. We dis- 
cussed, in turn, relative yields such as and K/h 
which allow insights into the hadronization conditions 
and help address questions related of chemical equilib- 
rium and non-equilibrium, as well as temperature of 
hadron freeze-out. 

Perhaps the most interesting result to pursue experi- 
mentally is the large value of 7s expected in hadroniza- 
tion of over-saturated QGP phase. The ratios, shown 
in Figs.|21 and ^ are not very sensitive to the choices 
we made that determine chemical potentials //b and fis^ 
they probe primarily the interplay between the 7s and 
T. We note further that most these ratios, at the fa- 
vored value 7s ~ 5, differ considerably from the chemi- 
cal equilibrium model expectations. The relatively large 
value of 7s we expect at LHC, twice as large as our 
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analysis finds at RHIC, derives from a larger absolute 
density of strangeness at hadronization of the decon- 
fined phase, combined with lower prevailing tempera- 
ture expected in deeper expansion supercooling. Indeed, 
while at RHIC 7?°^ I hadronization < 1, at LHC we expect 
7?°^ I hadronization ~^ 1.5-2. The difference in the phase 
space size of QGP with HG than leads to 7s — 5. 

We have further presented an in depth discussion of 
particle-antiparticle asymmetries which address the chal- 
lenge of chemical and strange quark chemical potential 



measurement. Since the particle-antiparticle yield differ- 
ence is small compared to each individual yield, a special 
effort will need to be made to acquire these difference 
ratios AN, ~ 0.1-4% at the level of 10% or better. 
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